Quantum feedback cooling of a single trapped ion in front of a mirror 
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We develop a theory of quantum feedback cooling of a single ion trapped in front of a mirror. By 
monitoring the motional sidebands of the light emitted into the mirror mode we infer the position 
of the ion, and act back with an appropriate force to cool the ion. We derive a feedback master 
equation along the lines of the quantum feedback theory developed by Wiseman and Milburn, which 
provides us with cooling times and final temperatures as a function of feedback gain and various 
system parameters. 
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I. INTRODUCTION 

Laser cooling and trapping of single ions 0,0,0 is onc 
of the highlights in the development of quantum optics 
during the last two decades. Single trapped ions are a 
laboratory paradigm of a quantum system, which can be 
prepared and controlled on the single quantum level, and 
whose time evolution can be monitored continuously by 
observing the scattered light in photodetection or homo- 
dyne measurements 0. By continuous observation of a 
single quantum system we learn the state of the sys- 
tem, as described by a conditional system density matrix 
p c (t), and based on this knowledge we can act back on 
the system, giving rise to quantum feedback control of the 
system of interest [EH EM El- I n the present paper we 
present a theory of quantum feedback cooling of a single 
trapped ion: by extracting from the scattered light the 
position of the ion in the trap, we implement a feedback 
loop on the system in the form of a damping force with 
the purpose of cooling the ion motion in the trap. Devel- 
opment of this theory is not only of fundamental interest 
in quantum optics, but the particular setup studied is 
motivated by ongoing experimental efforts |lllll2| to im- 
plement quantum feedback cooling of single trapped ions 
in laboratory. Indeed the present theoretical results pro- 
vides a quantitative basis for the understanding of these 
experiments |l2j. 

The particular setup studied in the present paper is a 
single laser cooled trapped ion in front of a mirror [T^. 
as illustrated in Fig. 2] and motivated by present exper- 
iments 0, ^2 ■ A single ion is stored a distance L from 
a mirror in a harmonic trapping potential. The ion is 
assumed to be a two-level system weakly driven near res- 
onance by laser light. Light is scattered into both the 
mirror mode, as well as the other other "background" 
modes of the radiation field. By detuning the laser on 
the red side of the atomic transition, the ion is laser 
cooled to a temperature corresponding to Doppler limit, 
where the mean occupation of the trap levels is much 
larger than one (i.e. far away from the sideband cool- 
ing limit to the ground state of the trap). Motion of the 
ion adds sidebands of the light scattered into the mir- 
ror mode displaced by the trap frequency. Observing the 



scattered light of these motional sidebands allows us to 
infer the position of the ion in the sense of continuous 
measurement theory, and feed back a damping force pro- 
portional to the momentum to implement quantum feed- 
back cooling. In this paper we will first formulate a con- 
tinuous measurement theory to read the position of the 
trapped ion from the scattered light using the language 
of stochastic Schrodinger Equations 0, U- Building on 
general quantum feedback theory formulated by Wise- 
man and Milburn [y, |7j , we will then derive a quantum 
feedback master equation for the motion of the trapped 
ion. This will allow us to study the dynamics and limits 
of quantum feedback cooling. 

For the setup studied in this paper the continuous 
readout of the ion position is based on light scattering 
into the mirror mode, with additional photons scattered 
into all other "background modes" of the radiation field. 
Spontaneous emission is intrinsically associated with a 
momentum recoil of the ion, which perturbs the ion mo- 
tion, i.e. contributes a heating mechanism for the ion. 
In a parallel paper |14| we study a quantum feedback 
scheme based on a dispersive readout of the velocity of 
the trapped ion to avoid this unwanted heating. It is 
based on the large variation of the index of refraction 
with the Doppler effect near a dark state resonance in an 
atomic A-system (based on electromagnetically induced 
transparency) . 

The paper is structured as follows. Sec. II presents the 
basic dynamic equations for the motion of an ion in front 
of a mirror. Quantum feedback equations are formulated 
in Sec. Ill, while results are presented in Sec. IV. 



II. MODEL AND BASIC EQUATIONS 

In this section we will develop the basic equations 
for continuous measurement of the photons in the mir- 
ror mode of the electromagnetic field. We will start 
with a detailed description of our model in terms of a 
Schrodinger equation for the coupled atom-bath system 
and the exciting laser. Continuous measurement theory 
provides us with a quantum stochastic Schrodinger equa- 
tion and hence a quantum stochastic master equation 
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in the Lamb-Dicke limit, where we adiabatically elimi- 
nate the excited state from the two-level atom. We will 
then derive the photocurrent obtained by detecting mir- 
ror mode photons and the corresponding stochastic mas- 
ter equation for the conditional density operator in the 
white noise (diffusive) limit. 



A. Single trapped ion in front of a mirror 

We consider a single trapped ion which is placed 
at a distance L from a mirror as indicated in Fig. ^ 
IH US- EE- 113- For the harmonic motion we assume 
a ID model in the z-direction (identical to the mirror 
axis). The harmonic trap has an oscillation frequency 
vt, and we denote the destruction (creation) trap oper- 
ator by a (a T ). The electronic degrees of the ion form 
a two-level atom with atomic transition frequency uj eg , 
with ground state \g) and excited state |e). We drive the 
two-level system with a laser with frequency lot, which 
couples the ground to the excited state with the Rabi 
frequency Q and a detuning from the atomic resonance 
= u>l — <-u eg . The atomic system Hamiltonian can 
thus be written as 
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H sys = v T a'a- A L |e) 



|e> (<?|+h.c. 



(2.1) 

Note that in this paper we set h = 1. In the interaction 
term we allow for a laser field incident at an angle x with 
respect to an axis normal to the z-axis. The momentum 
recoil due to absorption of a laser photon is represented 
by k e fiZ = 77 sin x (a + = f\ (a + a^) where the Lamb- 
Dicke parameter r\ = 27rao /A is the ratio of the size of 
the ground state and the laser wavelength. Due to the 
geometry of the system in consideration, the (quantized) 

electric field consists of two contributions, E^ = Em^ + 
E^' , where the Em denotes the modes restricted by the 

boundary condition of the mirror and E[ +) the remaining 
background modes [HG3. We adopt a ID model for the 
mirror mode and write for the electric field operator 



ElV (*) = i/ 
Jo 



dio a^e sin (fc(w)z) b m {io) (2-2) 



with a u a normalization factor for the mode function. 
The internal states of the atom couple to the vacuum field 
by an electric dipole transition. Denoting by d the dipolc 
matrix element, and introducing Pauli operator notation 
for the two level system, er_ = \g) (e\, the system-bath 
coupling Hamiltonian is 



'(f)(i)+4 + )(l)) ff _+hx. (2.3) 



The total Hamiltonian for the ion coupled to the radi- 
ation field is 




H — H K 



bath 



Hu 



(2.4) 



1 Laser 



FIG. 1: Physical Setup: L is the distance between the ion trap 
center and the mirror, \ is the incident laser angle. The light 
is collected in the photodetector PD. The feedback circuit 
consists of a bandpass filter, a phase shifter and an amplifier. 
The current Jfb is fed back to electrodes creating an additional 
potential for the ion. The mirror axis is equal to the z axis. 



Here -Hbath is the free Hamiltonian for the radiation field. 
We write this Hamiltonian as the sum of a Hamilto- 
nian for the mirror and the background modes -Hbath = 
H m + Hi,. In our ID model the mirror mode Hamilto- 
nian has the form H m — J dio Lob^ioj^m^) with b m (oj) 
photon destruction operators, satisfying bosonic commu- 
tation relations [b m (ui), b m (oj')] = 5(u> — lo'). Similar ex- 
pression can be given for the background modes. 

In analyzing this problem we are interested in the sit- 
uation where the time delay tm = 2L/c of the emit- 
ted light bouncing from mirror back to the atom is 
much shorter than the system time scales, in particular 
the spontaneous emission time from the excited state, 
tm <C 1/r, and the timescales associated with the laser 
interactions r 1/^2, 1/|Al|. This justifies the Markov 
approximation for the emission into the mirror modes, 
where we refer to 0] for a complete analysis. 

In the following we will denote the total spontaneous 
emission rate of the atom by T = T m + r^. Here Y m = eT 
with e the fraction of the solid angle covered by the lens 
is the emission rate into mirror mode, and = (1 — e)T 
the emission rate into the background modes. 



B. Quantum Stochastic Schrodinger Equation 

The dynamics of our model is summarized in the 
Schrodinger Equation 

|*(t)) =[-Lff sys + (2.5) 

+ Vf^cr- sin {k eg {L + z)) bljt) + h.c.+ 

, r+i , 

+ J du ^N{u)a^e-' mk ^ z bi{t) + h.c] |*(t)) 

We choose to formulate the problem in the language of 
a Quantum Stochastic Schrodinger equation (QSSE) 0, 
which allows for a direct connection with continuous mea- 
surement of the scattered light, and provides a direct link 
to quantum feedback theory developed in the following 
subsections. 
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In Eq. (|2.5() \^(t)) is the Schrodinger state vector of the 
combined atom-field system, i.e. the laser-driven trapped 
ion including the mirror and background modes of the 
radiation field. The first term on the RHS is the time 
evolution due to the system Hamiltonian (|2.1I) . 

The second and third line describe the interaction of 
the two-level atom with the mirror mode and the back- 
ground modes, respectively. We assume that these radi- 
ation modes are initially in the vacuum state. In writing 
Eq. (|2.5ll we have transformed to an interaction picture 
with respect to the free Hamiltonian of the radiation field 
Hb 0- As a result, we have introduced bath operators 
for the mirror mode b m (t) = l/V^r f dujb m (uj)e~ lut . In 
the Markovian limit these operators satisfy bosonic com- 
mutation relations 

[b m (t),bl(s)]=6(t-s). (2.6) 

In a Quantum Langevin formulation |5j b m (t) represents 
a quantum noise operator. Thus the second line of (|2.5|) 
describes the emission of photons by the atom into the 
mirror mode, with the center of the ion trap displaced a 
distance L from the mirror. We note that the motion of 
the ion couples to the light via the recoil, as seen by the 
appearance of z in the mirror mode function. This cou- 
pling imparts information of the ion motion on the light 
emitted in the mirror mode. In the following subsections 
we will analyze this scattered light to continuously mon- 
itor the atomic motion, with the goal of implementing a 
feedback loop to cool the ion. The coupling strength to 
the mirror mode is proportional to the square root of the 
spontaneous emission probability into the mirror mode 
r m = eT with e the fraction of the solid angle (typically 
e is much smaller than one). 

The third line in Eq. (|2.5|l represents spontaneous 
emission of the ion into the background modes. This 
is a coupling term familiar from the theory of laser cool- 
ing of two- level atoms 15j. Spontaneous emission into 
the background mode is again associated with a recoil of 
the ion motion. In our ID model for the motion of the 
trapped ion, it is the projection of this momentum on 
the z-axis which is the relevant momentum transfer. De- 
noting by 8 the angle between the emitted photon and 
the z-axis, and u = cos 9, we associate the transition 
for the excited state to the ground state including the 
momentum transfer with the operator e mfec92 (T_, where 
k eg = uj eg /c pa k^. Spontaneous photons can be emit- 
ted in all directions into the background modes consis- 
tent with the dipole radiation pattern of the given elec- 
tronic transition. We denote this (normalized) angular 
dependence by N(u). Thus the integral over u in the 
last line of Eq. 1(2.5(1 realizes photon emission into all of 
these possible directions. The operators b u (t) are again 
photon destruction (or noise) operators associated with 
these emission directions. They satisfy commutation re- 
lations 

[b u (t),bl(s)]=S(u-u')S(t-s), (2.7) 



and commute with the mirror bath operators b m (t) intro- 
duced above. The coupling stren gth to th e background 
modes is proportional to ^/Tb = y/O- — e)T. For red laser 
detuning < the cycle of laser excitation followed by 
spontaneous emission into the background mode leads to 
laser cooling. 

C. Ito form of the Quantum Stochastic 
Schrodinger Equation 

To give a meaning to the white noise limit (compare 
Eqs. I(2.6I2.7|) ). we must interpret the Schrodinger equa- 
tion (|2.5|l as a quantum stochastic Stratonovich equation 
[3 • As usual, it is more convenient to work with an Ito 
form, where Wiener noise increments 

/t+dt 
dsb m , u (s) (2.8) 

"point to the future" , i.e. are statistically independent 
of |\&(f)). These Wiener noise increments satisfy the Ito 
table 

dBmi^dB^t) = dt, (2.9a) 
dB u (t)dBl(t) = 5(u - u')dt, (2.9b) 

which follow from Eqs. (|2. 612.7(1 . the other entries of the 
Ito table being zero. The resulting Ito QSSE is 

(I) d|*(t)) =[-\H cS dt + ^TZC m (z)dBl (t) (2.10) 

, r+l , 

+ VnJ i du^/Wtu)C u (z)dBt(t)}\y(t)) 

Here, we have introduced the "jump operators" 

C7 M (£)=e- iufc -V_, (2.11a) 
C m (z) =sin(k eg (L + z))a-, (2.11b) 

which are associated with the emission of a photon in the 
background modes and the mirror modes, respectively. 
Furthermore, we have defined an effective non-hermitian 
system Hamiltonian 

H cS = H sys ~ i [r 6 + T m sin 2 (fc e3 (i + £))] |e) (e| . 

(2.12) 

The non-hermitian part of H c s arises from the Ito cor- 
rection in the conversion process. Physically, it corre- 
sponds to the radiation damping of the excited state due 
to the total radiation field. We also note that the photon 
absorption terms have disappeared in Eq. ((2.10(1 due to 
dB mjU {t)\^ (t)) = 0. This follows from our assumption of 
an initial vacuum state. 



D. Quantum Stochastic Master Equation 

We are interested in the time evolution of our sys- 
tem where the photons emitted in the mirror mode are 
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detected by a photon counter, while the background 
modes remain unobserved. Therefore, we are only in- 
terested in the dynamics of the reduced density opera- 
tor W (t) = Tr h {|* (*)) (* (t)\} where we trace over the 
background modes of the radiation field. We emphasize 
that W (t) still contains all the degrees of freedom of the 
mirror modes, in addition to the internal and external 
atomic dynamics. 

Using Ito calculus (see Appendix we obtain the 
quantum stochastic master equation (QSME) 



(I) dW (t) = - i [H cS W (t) - W (t) H f cS j dt (2.13) 
+T m J [C m (z)\ dBl (t) W (t) dB m (t) 

(c m (z)dBlW(t) + W(t)Cl(z)dB, 

+T b / duN(u)J[C u (z)]W(t)dt 



with H c ff defined in Eq. l|2.12[l . For the "recycling terms" 
we use the notation 



J [c] p — cpc^ . 



(2.14) 



Before proceeding we note that for e = 0, i.e. no cou- 
pling to the mirror modes, Eq. (|2.13|) reduces to the stan- 
dard master equation for ID laser cooling of a two-level 
atom pjj. In this case W is only an atomic density oper- 
ator containing the internal and motional dynamics. For 
e ^ 0, we still have a stochastic equation with the mirror 
bath degrees of freedom included. 



E. Adiabatic elimination of the excited state and 
Lamb-Dicke limit 

We will simplify the above QSSE fl2~TU|) and QSME 
(|2.13|1 with two assumptions. First, we assume weak 
laser excitation to the excited state, ft <C max(r,|A|). 
Second, we assume a small Lamb-Dicke parameter 7/ = 
2nao/X <C 1 (tight trap): this allows us to expand the 
exponents e [ki = e ir,< - a+a ^ = l + ir](a + a^ + 0{i] 2 ). Both 
of these assumptions are well satisfied in present experi- 
ments 

To eliminate the weakly populated excited level, we go 
back to Eq. (|2.10() and expand the state vector \^(t)) 
into ground state and excited state components, 



|* (t)) EE \i> g (t)) ® \g) + |Ve (*)) ® |e) 



(2.15) 



As shown in Appendix iBl we can eliminate \ip e (t)} in per- 
turbation theory in the Ito QSSE (|2.1()|l to obtain an 
effective equation for \ip g (t)). In a similar way as for 
Eq. (|2.13[) we obtain a QSME for the partially reduced 
density operator 



given by 
(I) dw(t) 



h eS w(t) - w(t)hl s 



dt 



(2.17) 



)] dBl(t)w(t)dB m {t) 
+ V7 (cm (z)dBl (t)w(t) + w (t)4 (z)dB m (*)) 
+ C b w(t)dt. 

The first three lines give the dynamics of the ion motion 
coupled to the mirror mode. The fourth line describes 
the traced-out action of the background mode on the ion 
motion, i.e. laser cooling of the ion. 

In Eq. H2.17fl we have defined an effective Hamiltonian 
acting only on the motional states of the ion, 

Kb = H t - ~74,(*)cm(*)- ( 2 - 18 ) 

where we expand the eliminated jump operators to sec- 
ond order in the Lamb-Dicke limit with the center of the 
trap at k eg L — 7r/4: 



Crn(z) « ( 1 + 1) (tt + (J) - (a + <J)^ 



The parameter 



7 = er- 



4 A| + E 



(2.19) 



(2.20) 



is the optical pumping rate into the mirror mode. The 
first three lines of Eq. 12.17fl thus describe the motional 
state coupled via laser excitation followed by spontaneous 
emission to the mirror mode. 

The Liouvillian L b in the fourth line of Eq. (|2.17|) 
is the standard laser cooling Liouvillian for weak field 
excitation and in the Lamb-Dicke limit [TL f^. fl5| ■ 

C b w (t) = A-V [a] w (t) + A+V [a)] w (t) (2.21) 
ee T cS {N + 1)V [a] w (t) + T cS NV [a f ] w (t) , 



where we have used the notation 

1 



T> [c] p ee epe) — — (crop + pc^c) 



(2.22) 



The rates 

A± = T) 



2 fl 2 
2 T h 



sin 2 x 



r 2 



(2.23) 

have the meaning of cooling (heating) terms for red laser 
detuning A L < 0. With r off = A_ - A + > and for 
Al < we have 



N 



A, 



A- - A. 



(2.24) 



w(t) = Tr b {\Tp g (t)) (t)\}, 



(2.16) 



which is the final mean trap occupation established by 
laser cooling (alone) . We have also used the abbreviation 
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e, n + 1 ) 




,n+ 1) 



n+ 1> 




\g,n- 1) 



FIG. 2: Contributing processes in the master equation at low 
intensity, (a) shows the cooling and heating terms due to the 
coupling of the laser to the motion and (b) shows the diffusion 
terms due to spontaneous emission [l8|. 



a = J du u 2 N(u) for the dipole transition parameter and 
X is the incident angle of the laser beam. With these 
definitions the mirror mode optical pumping rate (|2.20(l 
can be written as 7 = eNT c g/(l + a)rj 2 , and from F e g oc 
sin 2 x an d N oc 1/ sin 2 x we see that this pumping rate is 
independent from the angle of the incoming laser beam. 

In the following we will study a scenario where 
the laser cooling establishes a steady state with a mean 
trap occupation N 3> 1 (i.e. far from the ground state), 
as represented by the second line in Eq. I|2.17[l . This is 
the limit of Doppler cooling, which is obtained if T ^ 
vt- The minimally obtainable steady state energy in 
this limit is HT(a + l)/2. By observing the spontaneous 
emission into the mirror mode (see first two lines of Eq. 
()2.17(l ). we will infer the position of the atom to apply 
a feedback loop to cool the system (far) below the laser 
cooling limit. 

For completeness we note that in the case where the 
mirror mode is not observed, the reduced system density 
operator p(t) = Tr m {w(t)} obeys the master equation 



p{t) = -mT,p{t)] 

= -i[H T)P (t)} 



7 2? [c m (z)]p(t)+£ b p(t) 
£hcp(t) = C p(t), 



(2.25) 



which contains the dynamics from the free ion motion, 
and the dissipative dynamics from the emission into the 
mirror mode and laser cooling. In a second order expan- 
sion in terms of 77, we have 



V\ 



(z)] = rj 2 cos 2 (k eg L) V[a + o)]+0 (rf) (2.26) 



which, multiplied by 7, is typically much smaller than 
r e ff TV and thus the corrections in the heating and cooling 
rates will be neglected here. 



given by the internal and external degrees of the ion, for 
a given count trajectory is represented by a conditional 
density matrix p c (t) p| . 

Given the state of the system at time t, p c (t), the 
detection of a mirror mode photon in a time interval 
(t, t + dt] is associated with a quantum jump of the atom 
described by 



Pcjump (t + dt) 



J [c m ] Pc (t) 

Tr {J [c m ]p c (t)} 



(2.27) 



where according to l|2.11b|l the atom returns to from 
the excited state to the ground state, and momentum 
is transferred to the ion motion in accordance with the 
mirror mode function. In the case of no observed photon, 
the system evolves with the effective non trace-preserving 
Liouvillian Lq 



J c,no jump 



(t + dt) = (1 + L dt)p c (t) 



(2.28) 



where 



Lqp = -i 



h cS p - ph\ s 



CbP 



and ft- c ff is defined in Eq. I|2.18|l . The expected num- 
ber of counts in the interval (t,t + dt] is with dN c (t) = 
N c (t + dt) - N c (t) 

(dN c (t)) = p&*+*L = 7 Tr sys {J [c m ] p c (t)} dt (2.29) 

In view of dN c (t) — or 1, for this point process we have 
the Ito table dN 2 (t) = dN c (t) and dN c (t) dt = 0. 

We can summarize the above a posteriori time evolu- 
tion in an Ito stochastic Schrodinger equation (see. eg. 

0) 



(I) Pc (t)=C Pc(t)dt+ 

J [C m ] Pc (t) 



(2.30) 



Tr sys {J [c m ] p c (t)} 



Pc (t) 



x (dN c (t) - 7 Tr sys {J [c m ] p c (t)} dt) 

where £0 is defined in 12.25fl . This equation gives the 
time evolution of the conditional density matrix of the ion 
p c (t) for a particular count trajectory. Not observing, i.e. 
tracing over the mirror mode, is equivalent to taking the 
ensemble average over all count trajectories in (|2.30() . In 
this case, we recover the master equation p(t) = £op(t) 
for the a priori dynamics |5( . 



G. Diffusion approximation 



F. Continuous observation of the mirror mode 

We measure the photons emitted into the mirror modes 
by a photon counter as shown in Fig. ^ We denote by 
N c (t) the number of photon counts at time t. A particu- 
lar count trajectory is characterized the photon detection 
times t\, t2, ■ ■ ■■ Our knowledge of the state of the system, 



In the previous subsection we considered photon count- 
ing of the light emitted in the mirror modes, and the as- 
sociated time evolution of the system described by the 
condition density operator p c (t) . We are interested in 
learning the motion (position) of the atom from the scat- 
tered light in the sense of continuous measurement. The 
goal is to use this information to control the motion of 
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the atom, and eventually act back on the atom to cool 
it. 

The scattered light of a weakly driven trapped atom 
|l8j consists of (i) a strong elastic component at the fre- 
quency of the driving laser (see vertical transitions in 
Fig- EJ), and (ii) weak motional sidebands at the trap fre- 
quency vt suppressed by the Lamb-Dicke parameter rj. 
The information on the motion of the atom is encoded in 
the "motional sidebands" . We find it convenient to for- 
mulate the problem in a way, where we focus directly on 
the contributions of these sidebands to the photon count 
signal. The physical picture is that the elastic component 
acts like a "(strong) local oscillator" which beats with the 
"(weak) light emitted from the sidebands". This situa- 
tion is reminiscent of homodyne measurements in quan- 
tum optics 0, 0| , and will lead in the following to a de- 
scription in terms of a diffusive stochastic process rather 
than a point process associated with the photon count- 
ing described above. The formal expansion parameter is 
77 <C 1 (Lamb-Dicke limit). 

From the previous subsection we know that the mean 
number of photon counts in (t, t + dt] is 

(dN c (t)) kL = /4 lydt + 77? (a + a t > c (t)dt + 0{n 2 ). 

(2.31) 

The first term is elastic scattering. The second term, 
which is first order in 77, is proportional to z = a + , 
i.e. includes information on the ion motion. Here and in 
the following we take the center of the trap to be on the 
slope of the standing wave, i.e. k eg L = 

Following the analysis of homodyne detection 0,0, we 
split the stochastic variable dN c (t) into a deterministic 
and a remaining stochastic part, thus defining dY c (t), 



dN c (t) = pdt + T]dY c (i) 



(2.32) 



and we can show (cf. Appendix IcT) that dY c (t) is a Gaus- 
sian stochastic variable with non-zero mean, i.e. 

dY c {t) = dW(t) + 7 (S) c (t)dt 

with dW(f) a Wiener increment satisfying dW 2 (t) = dt. 

This leads us to define a photocurrent where we sub- 
tract the large constant contribution from the elastic 
scattering process, 



Ic(t) = T) 



dY c (t) 



dt 

IV («) c (t) + 



(2.33) 



In the diffusive approximation the conditional density 
matrix p c (t) Q, |j| obeys 



(I) dp&) 



Codt 



dW(t) H r , 



Pc(t) (2.34) 



where 



2f(t). 



H m Pc{t) = 277 (zp c (t) + p c (t)z -2(z) c (t) Pc (t)) (2.35) 

and Eq. (|2.34l) is derived from (|2.30() in the diffusive 
limit 77 -c 1 (cf. Appendix . 



III. QUANTUM FEEDBACK COOLING 

In the previous section we have reformulated the con- 
tinuous observation of the ion motion through sponta- 
neous light scattering into mirror modes in a form remi- 
niscent of homodyne detection. This will allow us below 
to study feedback cooling of trapped ions building on the 
Wiseman-Milburn theory of quantum feedback 

In Eq. 1|2.33[) we have obtained a current which is pro- 
portional to the mean value of the position of the atom. 
We want to use this information to feed back an appro- 
priate force proportional to the momentum to damp the 
motional state of the atom 0, ^^|. The information 
about the position is encoded in the motional sidebands 
of the current. In a harmonic trap of known frequency 
any combination of the average position and momentum 
can be obtained by shifting the sideband current by a 
phase of </3, if the trap frequency is much faster than any 
other (cooling) timescale in the problem (weak coupling 
limit). This phase (f> can be controlled electronically, and 
for cf> = 7r/2, the shifted current follows the momentum. 
A force, which is proportional to this current, can damp 
the motion of the ion. 



A. Feedback current 

We model the feedback circuit as shown in Fig. [3] 
First, the signal I c (t) given by Eq. (|2.33|l is mixed with a 
local oscillator of frequency uj ss vt to shift the signal of 
the motional sideband to zero frequency. Then the cur- 
rent is sent through a band pass filter of width B to cut 
off rapidly oscillating terms. The filter is described by a 
filter function Z(u>), centered around zero frequency. At 
the end the signal is mixed again with the local oscillator 
and amplified by a factor G. The feedback current can 
then be written as 



with £ (t) Gaussian white noise (£ (t) £ (£')) = 6 (t — t') 
(shot noise). We see that I c (i) follows (S) c (t) and thus 
represents a continuous measurement of the position of 
the ion. The information on the motion is contained 
in the sidebands of the current, i.e. in the frequency 
components centered around ±vt- 



Ifb,c (t) = Gcos(woi) / drZ(t — t) cos (lu t + cf>) I c (t), 
J — 00 

(3.1) 

where Z (r) is the Fourier transform of the band pass 
function Z(lo). The feedback Hamiltonian is specified in 
the next subsection. 
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FIG. 3: Electronic Setup of the Feedback loop as modeled in 
this paper. PD is the photodetector, LO is the local oscilla- 
tor, the B-box is the band pass (BP) filter and the 0-box is 
the phase shift. The LO signal is mixed to the current and 
subtracted after the BP filter. 



To evaluate the expression for the current, it is con- 
venient to change to a basis which is rotating with the 
frequency of the local oscillator loq by applying the uni- 
tary transformation U = cxp(— iujoa^at). The evolu- 
tion timescale of the density operator in this new frame, 
p c {t) = Up e (t)U' is determined by the detuning S = 
L0{) — vt and the cooling rates G"f,T e g. Under the as- 
sumption, that these frequencies arc smaller than the fil- 
ter bandwidth B, the feedback current is given by 



hh. 



(t) = G 



cos(woi)- (3-2) 



The first term in this expression, (X^) 1 = Tr sys {X0/5 c (i)} 
is the slowly varying expectation value of the quadrature 
component 



ae^ + a'e" 



(3.3) 



(in the rotating frame). The second contribution in 
Eq. I|3.2|l is defined as 



S(t) = dr cos {u) Q t + (j)) Z(t - r)£ (r) 



(3.4) 



It describes the noise which passes through the feedback 
circuit. The stochastic mean of E(t) is zero due to the 
vanishing mean of the white noise variable £(i), and the 
correlation function is given by 



(E(i)S(t')) « S B (t-t') + 0(- 



(3.5) 



Here <5b(^ — t') denotes a delta- function for functions 
which vary on a slow timescale much larger than B^ 1 . 
Thus for a clear separation of timescales, 



G7, S, T e ff < B < cj , v T , 



the current given in Eq. 113. 2[l is proportional to the slowly 
varying expectation value of Xq, and has a noise term 
which is delta-correlated on a timescale of the system 
evolution in the rotating frame. 



B. Quantum Feedback Dynamics 

The feedback current f (|3.2|l ) for <\> = —ir/2 is propor- 
tional to the slowly varying momentum of the particle. 
For the cooling of the ion motion, we apply a linear force 
which is proportional to the the feedback current (|3.2|l . 
For a trapped ion, this can be realized by applying a volt- 
age on the trap electrodes, which leads to a displacement 
of the trap center. The effect of the feedback force is 
given by the interaction picture Hamiltonian 



Hib = ifb, c (i - T)z/(i). 



(3.7) 



In this equation, zj(t) = U'zU is proportional to the 
position operator in the interaction picture, while r de- 
notes the finite time delay in the feedback loop, which 
we require to be much smaller than the trap frequency 
t -C 1/vt. The master equation i|2.34[l has to be com- 
plemented with the feedback term, 

(S) [dp c (t)]fb = - r) H) [2j(t), p c (t)] dt. (3.8) 

which has to be interpreted as a Stratonovich stochastic 
differential equation pj- For the slow dynamics of the 
density matrix in the rotating frame, we can make a ro- 
tating wave approximation and neglect rapidly rotating 
terms ~ exp(±2iw i)- The filtered noise l|3.4|l is delta- 
correlated on timescales slower than -B -1 , thus we have 
the coarse grained evolution of the density matrix 

(S) [dp c (t)U = ~ Tl (X^l (t - r)dtJC~ Pc {t) (3.9) 



G R 



dW E {t-T)fCp c (t), 



with the feedback operator 

K,p c (t) = -i [z,p c (t)] 



(3.10) 



and the "slow" Wiener increment dW-= (t) = S (t) dt. 

The total evolution of the system is determined by the 
conditioned master equation l|2.34|l plus the contribution 
from the feedback loop l|3.9(l . To combine the two equa- 
tions, we have to convert Eq. (|3.9|l from Stratonovich to 
Ito form. The total conditioned evolution is 



(I) dp c (t) = £ p c + x j ±HdW{t)p c {t) 



+ 7*7 (X^l (t ~ r)dt + ^ 7 Oi+ 
+ ^J\dW^t-T)\lCp c {t), (3.11) 



(3.6) where 



Cp c = ClcPc - ilSala, p~ c ] 



(3.12) 



(cf. Eq. ^2.'2b\i ) is the laser cooling Liouvillian in the 
rotating frame. 
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Because the exact photocurrent can not be kept track 
of in experiments, Eq. (|3.11|) is of limited use. The goal is 
to derive an equation for the ensemble averaged density 
operator. We follow the derivation given by Wiseman 
and Milburn in 7], where the measured current is fed 
back directly, and adopt it for our model. Assuming that 
the state at time t — t and all previous times is known, 
we take the ensemble average E[-\ of Eq. (|3.11|> over the 
trajectories in (t — T,t\. We then formally divide by dt 
and for convenience redefine p(t) = E[p c (t)]: 

(I) p(t)=£p(t) + |7»,<^)f(t-r)M*) (3-13) 

The density matrix p(t) is still conditioned on the evolu- 
tion up to time t — r, but not conditioned on trajectories 
in (t — T,i\. The ensemble average E[{Xi) c (t — r)p c (t)] 
factorizes because p c (t — r) is assumed known. Under the 
Markov approximation, we let r go to zero, while due to 
the coarse graining of the time evolution in Eq. (|3.8[) . dt 
will still be larger than this small delay An expansion in 
r yields 



PS) = [1 + O (r)] p c (t - t + dt) = 



[l + 0{r)\ 



1 



dW(t - t)H 



(3-14) 

Pc(t - T). 



We now can evaluate the remaining ensemble average in 
Eq. i|3.13[) because dW(t — t) is stochastically indepen- 
dent from p c (t — t). We obtain 

E [E(t - r)p c (t)} =VinE [H(f - r)£(t - r)} p(t) (3.15) 



-2{Xi) x e (t-T)p(t)), 

and thus the term in the last line, a conditional expecta- 
tion value, cancels with the second term on the right hand 
side of Eq. I|3.13|l . In going from the first to the second 
line in Eq. I|3.15[) we have dropped terms ~ exp(±ia;o£)- 
With this last step, we can finally evaluate Eq. (|3.13|) 
and write down the quantum feedback master equation 
(compare for the motional degrees of freedom: 



G 



G 2 



Cp+-f m K(Xtp + P X 4 ) + — 1 K'p. (3.16) 



The first term on the right hand side C is the laser cool- 
ing Liouvillian i|3.12[l in the rotating frame. The second 
term with K, given in Eq. I|3.10|l in the master equation 
is the feedback term. It acts back on the system and is 
responsible for cooling if we choose the parameters 5 and 
4> appropriately. The last term in the master equation is 
a diffusive term of the form of a double commutator. 



IV. RESULTS 

In the last section we have shown that for a separation 
of timescales <5, T e s <B<w ,i/r we obtain an uncondi- 
tioned (non-selective) master equation for the motional 
density matrix in the rotating frame. By inserting the 
definitions of C and K, the master equation reads 



p = -i5[ata, p] + AJD[a] + A + T>[al]+ 

G G 2 
- i j7??[z, X^p + pX^\ - ~^l[z, [z, p}]. 



(4.1) 



We have used the previously introduced variables z — 
a+a 1 and A = ae^+aU-^. In the first line of Eq. |0~TJ) 
we recover the master equation for laser cooling, with 
the corresponding heating and cooling rates A± given 
in Eq. (|2.23(l . The second line describes the effect of the 
feedback loop, where 7 = £AT c ff/(l+a!)77 2 is the emission 
rate in the mirror mode and G is the gain parameter 
amplifying the feedback current. The first term in the 
second line depends on the phase shift 4> and as we will 
show below, leads to the expected damping for (f> — —n/2. 
The second term arises from the noise in the feedback 
current and leads to a momentum diffusion, i.e. heating. 

We will derive solutions of the feedback master equa- 
tion (|3.1(jl) , which is bilinear in the position and momen- 
tum z and p z ■ It is convenient to use a Wigner function 
representation p| of the den sity matrix. This gives rise 
to a Fokker-Planck equation [13 for the Wigner function 
W{z, p, t) with dimensionless position and momentum 
variables x\ = z — z\J tuvt /2 and X2 = p — p z / \/2misT, 



dW(z,p, t) 
dt 



d 



(4.2) 



1,3 



d 2 W(z,p,t) 
dxidxj 



The Kij are independent of the phase space variables 
and Dij is diagonal, thus Eq. 14.2fl describes an Ornstein- 
Uhlenbeck process with drift matrix 



K = — I 1 

2 V 2G777 cos tf> + 25 1 - 2Grff sin <j> 
and the diagonal terms of the diffusion matrix 

r cff 1 



2N + 1, 2N + 1 + -G z j 



(4.3) 



(4.4) 



Here we have introduced the dimensionless detuning S = 
S/T e g and decay rate 7 = 7/r c fj normalized with respect 
to the width of the sidebands. The Gaussian Wigner 
function is uniquely determined by it's first and second 
position and momentum moments, and we will use the 
notation 



r—s\ 



W 



dzdpz r p s W(z,p, t), 



(4.5) 
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FIG. 4: Number expectation value for variable gain and solid 
angle fraction for r\ = 0.1, <f> = — tt/2, S = and N — 15. 



which equals the symmetric expectation value of the cor- 
responding operators. The bilinearity of Eq. (|3.16l) with 
respect to position and momentum gives rise to a closed 
set of equations for the first and second moments of the 
Wigner function individually and are given in Appendix 

El 

We are interested in the motional energy of the ion, 
which is related to the expectation value of the phonon 
number by E = huxi^a^ a) + 1/2). The expectation value 
for the number operator can be read off from the second 
moments of the Wigner function: 



(afla) ee (n) = (z*) w + <P 2 ) 



' w 



(4.6) 



We will calculate this quantity for different choices of 
parameters in the following subsections. 



A. Cold damping 

In this subsection we show results for <f) = — 7r/2 and 
(5 = 0, i.e. the center of the band pass filter is set exactly 
to the trap frequency. As derived in Appendix \B\ the 
number expectation value for the steady state in this case 
is given by 



(n) s 



N + ±777 (gjV - 1) G + l^G 2 
1 + 2r/jG 



(4.7) 



Taking the gain G = yields (n) = N, i.e. if we 
do not use the feedback current to influence the ion, the 
steady state occupation will be the one for standard laser 
cooling. We see that the slope of the occupation number 
is negative at G = 0, i.e. 



d (n) ss /dG\ G=0 = -JV(2N + l)/2 < 0, 



(4.8) 



and for G — > oo it diverges (note that in our model G7 
has to be smaller than B). Thus our theory yields a non- 
vanishing optimal gain G m ; n for which the occupation 




FIG. 5: Number expectation value for variable gain and 
steady state occupation number for r\ = 0.1, <j> = — 7r/2, 5 = 
and e = 0.01. N is given by the values at G = of the 
different curves, from bottom to top, N = 5, 10, 15, 20, 30. 



number has a minimum smaller than N, 

^/l + 8(2N+l)r] 2 j- 1 



Gn 



2rfi 



(4.9) 



Inserting this into Eq. I|4.7|l yields an expression for the 
minimal occupation number: 



A(2N - 1)77 2 7 - 1 + WT+ 8(2N + 1)^7 

(n) = 

\ /nun 16r]2 ^ 

(4.10) 

With increasing solid angle e we collect more informa- 
tion about the motional state of the system and hence the 
minimum (n) ss is expected to decrease, which is shown 
in Fig. With increasing e the optimal gain is decreas- 
ing, because the feedback noise term is growing with G 2 
while the damping term is linear in G. 

We show in Fig. [S] the decrease in the steady state 
phonon number with the gain. The relative decrease is 
larger with a higher laser cooling steady state phonon 
number N . For lower N, the mirror decay rate 7 oc N 
is smaller and thus we get less information about the 
motional state of the atom, which limits the feedback 
cooling. 

We will now expand (ri) in the limit of large (N 3> 1) 
occupation numbers. For a series expansion of l|4.10|l the 
formal expansion parameter is Ny/e, thus an expansion 
in the (usually also small) e is only possible for very low 
N. We make an expansion for large TV in the opposite 
limit (Doppler limit), while the condition Ny/e ^S> 1 has 
to be satisfied. N can be tuned with e.g. with the laser 
detuning A^. Then the minimal occupation number ap- 
proximately reads 



("U = y + 4 



1 



1 



Ne 



(4.11) 



which implies that for a sufficiently large collection angle 
the minimal obtainable phonon number is above N/2 and 
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thus feedback cooling alone cannot give a steady state. 
The reduction in the energy of the ion with time is due 
to the reduction in (p 2 ) w , while (z 2 )^ is constant, as is 
shown in the time evolution in Fig. Thus the Wigner 
function for the steady state will not be rotationally in- 
variant, but "classically squeezed" in the momentum di- 
rection. 

A phase space picture can demonstrate the action of 
the feedback on the system state (see Fig. |^a)). By 
feeding back a linear force / to the ion, we effectively 
apply a unitary operator of the form 



U(t) ~ exp (-ifxt) . 



(4.12) 



This operator acts as a momentum kick on a state with 
a magnitude proportional to the momentum, which we 
have chosen by setting <p = — n/2. The points in the 
Wigner function will tend towards the x-axis, while the 
diffusion term will counteract the feedback term, leading 
to a steady state Wigner function. 

The difference in the position and momentum vari- 
ance can be quantified; we will give an expression for 
the amount of "squeezing" , i.e. the ratio between the 
two half-axis of the error-ellipse for the Wigner function 
in phase space is obtained by rotating the axes of the 
ellipse: 



semiminor axis 1 — / 
semimajor axis 1 + / 



f = 



' pp 



4cr 2 



(4.13a) 
(4.13b) 



Here a zz = (z 2 ) 



w 



' pp 



i w and (Jpp 



(P 2 )w-(P) 



are 



the variances of position and momentum, respectively, 



and a. 



( Z P)\ 



i w 



(P) 



w 



As mentioned, due to 



the affection of only the a pp component, a zp — in the 
case 4> = —tt/2. The range of the squeezing parameter is 
< r a < 1, where a small value corresponds to strong 
squeezing and for r a = 1 the state is symmetric. 

The time dependent Fokker-Planck equation is solvable 
analytically and the timescale of the cooling process is 
given by the eigenvalues of the drift matrix l|4.3|) , which 
are in this case T c ff and T e g + 2rffG corresponding to the 
usual Doppler cooling and the feedback cooling. This 
shows that the feedback cooling happens on a timescale 
faster than laser cooling alone. 



B. Variable feedback phase 

For a phase 4> ^ —tt/2, the magnitude of the feed- 
back force is proportional to the projection of the mo- 
mentum on an other rotated axis in phase space. We 
have pointed out in Eq. I|4.12l) that the action of the lin- 
ear force (shifted trap) is always a momentum kick. Thus 
the particle will always be "kicked too hard" or not hard 
enough towards the phase space center. We will calculate 
the regions of stability where the feedback can still lead 



FIG. 6: Time evolution of the variances for N = 15, 8 = 0, 
£ = 0.006, tj = 0.06, <t> = -tt/2 and G w 1.51 The solid line is 
the variance of the position (z 2 )w and the dotted line is the 
variance of the momentum (p 2 )w- 




FIG. 7: Number expectation value for different feedback 
phases. Here 8 = 0, thus the expectation values for (f> —* n — <f> 
yield the same results. Other parameters: N — 15, e = 0.01 
and rj = 0.1. 



to a steady state. Such a steady state will only form if the 
both eigenvalues of the matrix k are positive. One eigen- 
value of this matrix is r e ff for arbitrary <j>, giving again 
the usual Doppler cooling, and the other eigenvalue is 
r c ff — 2G7/7 sin 0, which is always positive for negative 
angles. For positive angles > 0, the gain has to ful- 
fill the condition G < T e g j '2 777 sin cf>. If this condition is 
satisfied, a steady state number expectation value exists 
and reads: 

( n ) ss = [(! - ?77Gsin0) (1 - 2777GSU1 ^p 1 x (4.14) 
x [N + - (4JV - 1)7776? sin </3+ 

+ -jG 2 (l + 4 7 ?7 2 (2N + 1 - 2 sin 2 <j>)) - 
8 

- ?77 2 G 3 sin(/)}. 

From Eq. (|4.14|) we can see that an energy decrease 
via feedback cooling is only possible for angles —tt < 
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12 3 



FIG. 8: Time evolution of variances {z 2 )w, {p 2 )w and (zp)w 
for N = 15, 8 = 5, e = 0.006, r) = 0.06, <j> « -85 degrees and 
G ~ 2.03 The solid line is the variance of the position (z 2 )w, 
the dotted line is the variance of the momentum {p 2 )w and 
the dashed line is (zp)w- 




FIG. 9: Schematic drawing of the Wigner function error el- 
lipses for the initial thermal state (circle) and the feedback- 
cooled state (ellipse). In (a), 4> — —tt/2 and 5 = 0, note 
that the position variance stays constant while the momen- 
tum variance is decreased. The action of the force is always a 
kick in the momentum direction and the force is proportional 
to the averaged momentum. In (b), 5^0 and <f> ~ —it/2, 
here the Wigner function is rotating and both variances are 
damped, resulting in lower energies. 



<j) < by calculating the slope d(n) ss /dG\c=o- Because 
Eq. H4.14fl is of higher order in G than the equation we 
had for <f> — —tt/2, l|4.7|l . we will not give an analyti- 
cal solution for the minimal gain and number occupation 
here. We also find that for (f> ^ —tt/2 the optimal occu- 
pation number is higher than <j) = —tt/2 (compare related 
studies in |2Cj ) . The steady state occupation number for 
varying cj> as a function of the gain is plotted in Fig. {7\ 
where we can see that for non-optimal phases the range 
of G for (n) ss < N is shrinking. 

For the special case of <j> = n or <j> = 0, no cooling 
can be observed any more and the number expectation 
value is quadratic in G. In principle, a steady state with 
(n) > N always exists with 

(n) ss =N+ IjG 2 (1 + 4r? 2 7 (2N + 1)) . (4.15) 
8 

The more interesting feature of the <p — tt case is that 
in the master equation (|3.16|l the feedback term (sec- 
ond term) reduces in a rotating wave approximation to 
a Hamiltonian term of the form — i lAva^a, p\ . For this 
case we observe a small shift Aj/ in the frequency of the 
trap linearly proportional to the gain. In this paper, we 
have not discussed the detailed experimental setup used 
to apply the force to the ion, which would be necessary 
for the knowledge of the exact forces acting on the ion. 
For (j> = 7r one can measure the frequency shift in the 
location of the sideband and determine the conversion 
factor from the gain parameter G used in this paper and 
an experimental gain factor, which might be the real elec- 
tronic gain in the feedback loop. 

C. Rotation in Phase Space 

We have shown that the phase <j> = —tt/2 we chose 
leads to the lowest energy of the motional state of the 



ion. The variance for the position operator (z 2 )^ = 
(2N + 1) /4 remains constant with time as shown e.g. in 
Fig. El thus posing a lower limit to the obtainable energy. 
The detuning 6 of the local oscillator in the feedback loop 
from the trap frequency creates a tunable slow rotation of 
the (interaction picture) Wigner function in phase space. 
This results in "squeezing" of all quadrature components 
(see Fig. Ett>)), and the Wigner function can regain a 
symmetric shape. Of course the timescale for this ro- 
tation has to be much slower than the filter bandwidth 
B. 

For the time evolution of the variances, the effect of 
the detuning is illustrated in Fig. |H1 We see the time 
evolution of an initially thermal (symmetric) state with 
an occupation number of TV. In contrast to Fig. the 
width of the Wigner function in the momentum and the 
position space are alternately decreased until they reach 
the new feedback steady value. For a larger detunings the 
two variances decrease equally in time and energetically 
lower states can be reached. For a rotation of the Wigner 
function with the frequency 5, we have to compare this 
rotation timescale with the cooling timescale 7. For 7 
comparable to (5 the optimal phase is is shifted with re- 
spect to —tt/2 because the Wigner function is rotating in 
phase space during the cooling time. When the detuning 
is much larger than the cooling rate, the Wigner function 
ellipse direction will not be resolved during the cooling 
time and thus the optimal phase returns to —tt/2. By nu- 
merical optimization (Fig. HOfl we find that the optimal 
phase is shifted from —tt/2 asymmetrically with respect 
to the detuning S. It reaches it's maximum excursion for 
a value of 5/r e g w 1, for higher detunings the optimal 
phase approaches —tt/2 again. For these optimal values, 
we plot in Fig. ^2 the squeezing parameter r a l|4.13a|l . 
which is one for a symmetric Gaussian state. We see 
that the state at no detuning is "classically squeezed" as 
we already mentioned in subsection llV Aj and the squeez- 
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and get for the occupation number: 

_ N-\rnG+\*,(? 
K /ss I + 7/7G 

The minimal occupation number for the same limit we 
took in deriving Eq. (|4.11|l we get for N ^> 1: 



(4.17) 



1 



2s 



1 _ 2(1 + a) 

2 8eN 



(4.18) 



This expression does not include the large term N/2 any 
more and thus the obtainable energy for large TV has an 
upper bound which is independent of TV, thus feedback 
cooling alone can give a thermal (symmetric) state with 
a temperature below the Doppler temperature. 



FIG. 10: Optimal phase as function of the detuning with 
77 = 0.1, e = 0.006 and N = 10 (solid line), TV = 17 (dotted 
line) and TV = 24 (dashed line). The weak dotted line marks 
4>= -it/2. 
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FIG. 11: Squeezing parameter ra for the optimal phase at a 
given detuning and the minimal energy state with e = 0.006, 
77 = 0.06 and with the variable parameter TV = 10 (solid line), 
TV = 17 (dotted line) and TV = 24 (dashed line). 



ing increases up to S « 1, then upon approaching 5 — > 00 
the squeezing parameter approaches one, and the state 
is thermal. 

For increasing 5, we also show that the number expec- 
tation value is decreasing. We will not give an analytic 
expression for (ri) for an arbitrary 5 here. We merely 
calculate the minimal number of phonons in the limit of 
large S. For this we require an additional separation of 
the timescales between the effective feedback cooling rate 
and the detuning, while the other timescale inequalities 
still hold: 

7<^<B. (4.16) 

With these new conditions we take the detuning S — > 00, 
where the optimal feedback phase is again <fi = —n/2, 



V. CONCLUSION 

In this paper we have studied quantum feedback cool- 
ing of a trapped ion in front of a mirror. This work is 
motivated by recent experiments [Tlj , and - as shown in 
|12| - provides a quantitative understanding of the ex- 
perimental results. 

In the setup discussed in this paper the final temper- 
atures are limited by the collection efficiency, s, and the 
constant scattering of photons for the position measure- 
ment. This combination of heating due to the recoil, 
and laser cooling due to the red detuning of the laser 
leads to a steady state temperature (Doppler limit). The 
effect of quantum feedback cooling is studied as an ad- 
ditional cooling mechanism on top of the ongoing laser 
cooling. For the experimentally relevant parameters this 
leads to sub-Doppler cooling, but it seems difficult to 
achieve ground state cooling in the t rap along these lines. 
As shown in a parallel publication [lj|, we can devise a 
purely dispersive and thus non-invasive readout of the 
velocity of the trapped ion based on the variation of the 
index of refraction with velocity, i.e. based on electro- 
magnetically induced transparency. Such a scheme al- 
lows, under idealized conditions, ground state cooling of 
the ion purely by quantum feedback. 
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APPENDIX A: DERIVATION OF THE 
QUANTUM STOCHA STIC MASTER 
EQUATION 

Starting from the stochastic master equation (|2.1U|) 
we define the reduced density matrix W(t) = 
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Tr b {|*(<)) (*(*)]}. Note that W(t) is now a trace-class 
operator for the internal electronic, the motional and the 
mirror mode bath degrees of freedom. We calculate 

dW(t) = Tr b {|*(f + dt)) (*(t + dt)\- |*(t)) (*(*)]} 

(Al) 

by inserting + dt)) = |*(t)> + d|*(t)) from 
Eq. Using the Ito rules dB u (t)dBl(t) = S(u - 

u')dt, and cyclic property of the trace for background 
bath operators, all terms of the form 

Tr b {dBt(t)m))m)\}= (A2) 
= Tr b {|*(t)> (9(t)\dB a (t)} = 

vanish because the initial bath state is the vacuum state. 
With these rules we obtain Eq. (|2.13() . 



APPENDIX B: ADIABATIC ELIMINATION, 
LAMB-DICKE LIMIT AND LASER COOLING 



with 

dC\ = \/r~ b J du y/N(u)dBl + x/TZsm(k eg L) dB^, 

(B3) 

dC\ = -iVTft J du y/N {u)udBl + cos{k eg L) dB^. 

(B4) 

Consistent with the above approximations we neglect 
here and in the following terms oscillating at twice the 
trap frequency vt- Physically speaking, the fourth line 
of Eq. I|B2|1 will correspond together with third line to a 
heating and cooling term, and the last line describes a 
diffusive term (cf. Fig. 

Taking the trace over the background modes to define 
a reduced density operator w(t) according to (|2.16ll we 
use the Ito rules, e.g. 

Tr b {dBl (t) |V S (t)) (ip g (t)| dB u , (t)}=S(u- «') p (t) dt 
to derive Eq. (|2T7|) . 



This appendix fills in the details of deriving the QSME 
(|2.17|) from the QSSE l|2.10|l under the assumption of 
weak driving and small Lamb-Dicke parameter. Note 
that we will need to consider two different Lamb-Dicke 
parameters due to the exciting laser which is not collinear 
with the z-axis. As in section Til Al we denote fj = r\ sin \. 
Inserting the ansatz (|2.15(l into the QSSE p.lOjl and 
transforming to an interaction picture with respect to 
Ht we get 



\i>e (*)> = 



if2 / 1 — ^fj 2 a^a 




i (A L - v T ) + t -i (A L + v T ) + 4. 



We insert this expression back into (|2.10() . We obtain 



d\i>g (*)) 



' 4 



-iA 



O e 



±1\U T t\ 



fj 2 a^a 



+ 



fj 2 aa' 



i (A L - vt) + 5 i (Al + v t ) + ^ 

in i 



(B2) 
dt+ 



-iA, 



lr/ae 



dC[- 



-lUTt 



ifja^e iVTt 



i(A L -v T )-j i{A L + v T )-j 



dC{- 



-iA L + § 



-dCl 



l^(*)> (I) 



APPENDIX C: HOMODYNE 
PHOTODETECTION AND THE DIFFUSION 
APPROXIMATION 



As we have seen in Sec. Ill Fl the statistics of the de- 
tected photons in the mirror mode are determined by the 
Poissonian stochastic variable dN c (t). Like in homodyne 
detection, where a strong local oscillator beats with the 
photodetection signal from a quantum system, an elastic 
scattering term beats with the signal given by the cou- 
pling of the light to the ion's motion (cf. Eq. i|2.31|) ). 
(Bl) The parameter which gives the difference in the magni- 
tudes of these terms is the Lamb-Dicke parameter rj. We 
split the stochastic variable dN c into a constant (deter- 
\ip g (t)) . ministic) part and a remaining stochastic part: 



dN c (t) = ^dt + r]dY c (t). 



(CI) 



The stochastic expectation value of this equation is al- 
ready known from Eq. (|2.31l) : 

(dY c (t))= 7 (z) c (t)dt. (C2) 

We check the distribution properties by calculating 

(dN c (t) - \^dt\ 2 _ dN c (t) 



dY c ' (t) 



\-ydt + T)dY c (t) 7 



(C3) 



which tells us that the stochastic variable has Gaussian 
properties, and thus is associated with a white noise prob- 
ability distribution. Thus dY c (t) = y/ r y/2/r]dW(i) where 
dW(t) is a Wiener increment. 



14 



The evolution of the system conditioned on measur- 
ing the photocurrent can be seen by expanding the first 
bracket in the stochastic master equation (|2.30(1 to first 
order in the Lamb-Dicke parameter rj, and noting that 
the second bracket in l|2.30[) is 



dN c (t) - (dJV c (i)> 



(C4) 



Thus, using the formal derivative £(t) = dW(t)/dt we 
obtain the conditioned equation for the reduced density 
matrix, Eq. I|2.34[) . 



APPENDIX D: EQUATIONS OF MOTION FOR 
THE MOMENTS OF THE WIGNER FUNCTION 



In Sec. IIVI we use a Wigner function representation 
for the density matrix and get an Fokker Planck equa- 
tion (|4.2[) equivalent to the master equation Eq. (|3.1t)|) 
with the drift matrix (|4.3|) and the diffusion term l|4.4|l . 
The equations of motion for the first and second mo- 
ments of the Wigner function in terms of the normalized 
position and momentum variables z = x\ and p = X2, 
respectively, are: 



d 



-q! {x k xi) w =2D ki + 2D ik - 



(Dl) 
(D2) 



[ K kj (%i x j) w + ^ (%k%j)^ 



For a constant drift matrix, the equations for the first 
moments are trivial, and if the eigenvalues of k are posi- 
tive, the steady state value is zero for both moments. We 
will therefore not concentrate on the first moments. We 
will give the equations for the second moments which are 
relevant for the number expectation value, and for this 
purpose we define a vector of second moments 

y(*) = ((z 2 (t)) w , (P 2 (t)) w , (mp{t)) w ) T ■ (D3) 

We can write the equation of motion in a compact form 
as 



y(t)=My(i)+u 
where the evolution matrix is 



(D4) 



M 

r e ff 



i o s 

l-2Gsin0 Gcosq^ + S 

2Gcos</> + 2<5 -25 1-Gsin0 



Here G = Grjj and 



(D5) 



(D6) 



The steady state results are obtained by setting y(t) = 0, 
which yields 



y ss = —M x u 



(D7) 



and we can calculate Eqs. (|4.7I) and (|4.14|l with (n) = 
Vi +V2- 1/2. 
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